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This paper describes a lattice version of the Skyrme model in 2 + 1 and 3 + 1 dimensions. 
The discrete model is derived from a consistent discretization of the radial continuum problem. 
Subsequently, the existence and stability of the skyrmion solutions existing on the lattice are 



m 



C/^ \ investigated. One consequence of the proposed models is that the corresponding discrete 

skyrmions have a high degree of stability, similar to their continuum counterparts. 
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Many attempts have been made in order to obtain the discrete analogue of given continuum 



t^^ ' systems that admit solitary wave solutions, so that their characteristics in the continuum 

0|^ . be preserved in the lattice. One such characteristic for the topological systems is often the 

^ ! existence of the so-called Bogomolnyi bound [1]; while for the integrable ones it may be 

rS ' the existence of the Lax pair [21. In the continuum, the stability of the topological solitons 

H : 

is often related to the existence of the energy bound while, the stability of the solitons in 
the integrable systems is related to the presence of an infinite number of conservation laws. 
However, in the lattice, only in a few cases the Bogomolnyi bound has been preserved [21 HI [5] 
while it is even more difficult to define the lattice Lax pair (especially in higher dimensions) of 
the corresponding integrable system [2J. Moreover, the stability of the solitons in question is 
not guaranteed. Often, under discretization, the topological properties are essentially lost (as 
may be expected since topology is related to continuity). Lattice versions of these nonlinear 
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wave bearing systems have been much studied (for purposes of numerical simulations [6J or 
regularization of the quantum field theory, or because they are of fundamental physical interest 
in their own right [7j). 

For a given continuum model, there are many different discrete analogues which reduce to 
it in the continuum limit. The object of this paper it to present a lattice version of the Skyrme 
model in 2 + 1 and 3 + 1 dimensions. The Skyrme model is a popular model of the dynamics of 
pions and nucleons, incorporating the former as its fundamental pseudo-Goldstone field and the 
latter as topological solitons. Its continuum version has been widely studied using numerical 
and analytical methods (for more details, see for example, Ref. |8]). Its lattice formulation 
is of some importance since the model is non-renormalizable in perturbation theory and thus, 
existing treatments of the model are semiclassical (quantizing only the collective degrees of 
freedom of the soliton). Full quantization of the theory requires a cutoff which can be attained 
by its lattice version. 

The most interesting feature of the Skyrme model is the stability conferred on the soliton 
by the topology [9]. An open question is whether this stability is preserved on the lattice. 
In this paper, lattice skyrmions obtained within an appropriate discretization of the Skyrme 
model in 2 + 1 and 3 + 1 dimensions are typically found to be stable in our parametric 
investigations. This suggests that these discretizations bear some important features of their 
continuum counterparts, while being easier to handle from a numerical point of view. 

Our presentation is arranged as follows. The next section reviews the baby Skyrme model 
and then reparametrizes the fields to impose radial symmetry. Only then is the model dis- 
cretized, as shown in section 2.1; while in section 2.2, numerical results on the the existence 
and stability of a single (radially symmetric) baby skyrmion on the lattice are given. The case 
of soliton configurations in the Skyrme model in 3 + 1 dimensions was dealt with in the same 
way in section 3. Finally, our conclusions are presented in section 4. 



2 The Baby Skyrme Model 

Let us begin with a brief review of the Skyrme model in 2 + 1-dimensions (so-called baby 
Skyrme model). The Lagrangian density of the model is of the form 

L = ^d^cP 9°0 - ^(9„0 X 9^0)(9"0 X d^cP) -^^l-{n- ct>f) . (2.1) 

The field is a map from M^ -^ S"^, where M^ is the 3-dimensional Minkowski space with 
metric 1]°'^ = diag(l, —1, —1) and the target space S"^ is the 2-sphere of unit radius embedded in 
Euclidean 3-space. Therefore, the field (/> is a scalar 3- vector with norm one, i.e. |0p = 1. The 
constants n, fi are free parameters which have the dimension of length and energy, respectively. 
The first term in (12.11) is the familiar 0(3) sigma model, the second term is the 2-dimensional 
analogue of the Skyrme term and the last term is the potential. 

The presence of the potential in (12. ip ensures the existence of stable skyrmion solutions. 
There are two other options in literature for its form: i) the holomorphic model where the 
potential term is (1 + n • (pY but stable skyrmions cannot be obtained and (ii) the old baby 
Skyrme model where the potential term is {1 — n ■ (p) and stable non-radially symmetric 
skyrmions exist. The Lagrangian (12. ip corresponds to the new baby Skyrme model [IQj which 
possess radially symmetric skyrmions. In [llj, approximate analytic skyrmion solutions of the 
new baby Skyrme model were obtained by exploring its topological properties. 

Finiteness of the energy requires the potential term to vanish at infinity, implying that 



linir^oo 0(i, 3^, y) = n (where r = yx^ + Z/^)- For simplicity we choose n to be the vacuum 
state, that is n = (0,0,1). Thus a topological number exists since the field 0, due to the 
boundary conditions, can be considered as a map from S"^ -^ S"^. 
The topological charge is the homotopy invariant of the field 

deg[(f>] = —/(/)• (d^cf) X dy(l>) dxdy. (2.2) 

and thus, conserved. 

2.1 Discrete Baby Skyrmions 

In order to obtain the discrete analogue of the baby Skyrme model we restrict our consider- 
ations to fields which are invariant under simultaneous rotations and reflections in space and 



target space. Thus, we assume that the field 4> is of the hedgehog form 

(l)i = kising{r,t)ki, (p^ = cos g{r,t) (2.3) 

where /c, for i = 1, 2 is a unit vector given in terms of the azimuthal angle 6 and the topo- 
logical charge N = deg[0] as ki = [cos N6, sin NO); and g{r,t) is the real profile function 
which satisfies certain boundary conditions. Then, the respective energy functionals related 
to Lagrangian (12. ip are 
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^kin = TT rgUl + —^ sin'g dr (2.4) 



f { 2 '^'^' 2.2 N^ ■ 2 2 . 2 A ^ 

-c^pot = TT [rg^-l g^sm g -\ sm g + firsm g j dr. (2.5) 

The boundary conditions for the skyrmions are: g{0,t) = tt and g{r,t) = as r ^ cxd. Note 
that the first and third term in (12.51) corresponds to the static 0(3) sigma model energy; the 
second term corresponds to the 2-dimensional static Skyrme energy and the last one is the 
potential. 

Hereafter, r becomes a discrete variable with lattice spacing h. So, the real-valued field 
g{r,t) depends on the continuum variable t and the discrete variable r = nh where n G Z^. 
Then, g^ = g{{n + l)h,t) denotes forward shift and thus, the forward difference is given by 
^9 = (fi'+ ~ 9)/h- There are many possibilities for discretizing the above energy functionals. 
However, based on the approach introduced in [5] for the discretization of the 0(3) sigma 
model, we assume that 

2f{h) . fg+-g 
sm 



h \ 2 

The parameter f{h) is an arbitrary function of the lattice spacing subject to the constraint 
f{h) ^ 1 as /i ^ 0. 

The origin must be treated in a special way since the functionals are not defined at n = 0. 
One possibility is to assume (following [5]) that at the origin we have 
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fco.= (^l. (2.7) 



Then, the kinetic and potential energy of the discrete baby Skyrme model are defined by 
the following expressions 
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Note that, for discretizing the second term of the energy at the origin, a combination of 
formula (12. 6p and (12.71) has been used. Actually, when this term is absent no stable lattice 
baby skyrmions can be obtained. Recall, that this term is the discrete analog of the Skyrme 
term which stabilizes the solution and apparently its presence is vital (even at the origin). 

For K = fj, = the model (12.81) becomes the discrete version of the 0(3) sigma model 
introduced in [5]. 

The lattice equations of motion obtained from the variation of the Lagrangian L = i?kin — 
Epot given by ([MD are 
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(2.9) 

Next, our task is to study whether the aforementioned lattice equations admit skyrmion solu- 
tions and if this is the case investigate whether the topological properties of the solutions are 
maintained in the presence of the lattice. 

2.2 Numerical Simulations 

Our numerical procedure for obtaining the baby skyrmions is the following: we use a fixed 
point iteration to identify the static solutions of equations (12. 9p . An initial guess (for the 
fixed-point Newton iteration) in the form of an inverse trigonometric function (an arccos in the 
radial direction) is used which subsequently, after a few iteration steps, converges to an exact 
stationary solution. Examples of such solutions are shown in Figure[T]for h = 0.65 (left panels) 
and h = 1.5 (right panels). The results are obtained for the choices: N = K, = ^ = f = l 
(unless noted otherwise), although variations of the parameters do not significantly affect our 
conclusions presented below. 

As expected, increase of h renders the solution more "coarse" (i.e., with a fewer sites 
participating in the "spine" of the skyrmionic structure). However, the stability of the wave is 
not crucially affected by the discretization, since the baby skyrmion remains stable throughout 
the range of parameters used in our numerical investigations. The stability is inferred from 
the eigenvalues of the relevant Jacobean of linearization around the solution. In particular, a 
linearization ansatz of the form: 

9n = Qn +eexp{\t)wn (2.10) 

leads to an eigenvalue problem (to 0(e), where e is a formal small parameter) for the eigenvalue- 
eigenvector pair (A,U7„); g^ represents the stationary solution obtained in the aforementioned 
Newton step. A solution is deemed to be stable if none of the eigenvalues A = A^ + iXi of 
the linearization problem is found to have a strictly positive real part A^. A particularly 
peculiar feature of the model is that the second term on the left hand side of equations (12. 9p 
does not contribute to the linearization (its lowest order contribution is 0(e^)). This suggests 
that even if the solution is found to be linearly stable, the full dynamics of the system (12. 9p 
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Figure 1: Typical profiles and stability results for the discrete baby skyrmions. The profile 
of the solution is shown in the top panels for h = 0.65 (left) and h = 1.5 (right). The 
corresponding spectral plane (Ar, Aj) of the eigenvalues A = A^ + ^Aj of the linearization around 
the solutions is shown in the bottom panels. 



should be considered, as it can, in principle, lead to nonlinear instabilities that cannot be 
detected at the linearization level. The results of the branch of solutions and their stability, 
obtained as a function of the lattice spacing h are summarized in Figure [2l The potential 
energy of the solutions, evaluated based on (12.81) . reveals that the discretizations contributes 
towards decreasing the energy of the stationary solutions. The eigenvalues vary as a function 
of h, and interesting features such as bifurcations of internal modes (see e.g. [121 [IS] and the 
references therein) arise e.g. for h > 1.15 from the bottom edge of the continuous spectrum. 
Additionally, a point spectrum (isolated) eigenvalue exists, which is also clearly discernible in 
the plots. Nevertheless, these do not appear to significantly affect the stability of the obtained 
discrete baby skyrmion structures. 

Finally, we examine the dynamical evolution of the baby skyrmions to confirm the stability 
of the solutions obtained herein. We initialize the full dynamics of (12. 9p with an exact baby 
skyrmion (for h = 0.65 in Figure [3]) perturbed by a random (uniformly distributed) field of 
amplitude 5 x 10"'^. In Figure [3] the space-time evolution of the waveform and its persistence 
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Figure 2: The left panel of the figure shows the potential energy of our exact static solutions 
given as a function of the lattice spacing h. In the right panel the principal (minimal, as well 
as maximal) imaginary parts of the eigenvalues of the linearization around the solutions are 
shown as a function of h. 



in time (left panel) is presented, and it is obvious that despite the strength of the perturbation 
(which is, however, substantial) remains bounded during its time evolution in the right panel 
of the figure. This indicates that the baby skyrmions are robustly stable dynamical structures 
of the corresponding discrete equations. 

Note that in [TJ], discrete 2-dimensional topological skyrmions have been constructed for 
a novel lattice version of the baby Skyrme model (i.e., the 2-dimensional topological Heisen- 
berg model) but their stability was far weaker than the continuum ones since a fairly small 
perturbation caused their decay. 



3 The Skyrme Model in 3 + 1 Dimensions 

The Lagrangian of the SU{N) Skyrme model in (3 + 1) dimensions can be written in terms 



of the currents R^ = d^UU ^ as 



Uii^C = -Ui{R,R'') - ^tr ([i?^,i?,][i?^i?'^]) 



(3.11) 



where we have used scaled units of energy and length, and a (+, — , — , — ) signature for the 
space-time metric. The asymptotic value of the SU{N) Skyrme field U{x,t) has to tend to a 
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Figure 3: Space-time evolution of the baby skyrmion for h = 0.65 under an initial (random, 
uniformly distributed) perturbation of amplitude 5 x 10^'^. Note that despite the significant 
strength of the perturbation, the solution profile remains essentially intact (left panel), while 
the difference 6g{n, t) = g{n, t) — g{n, 0) remains bounded (right panel) attesting to the robust 
stability of the discrete baby skyrmion. 



constant matrix at spatial infinity, so that finite-energy configurations can exist. 

The baryon density, whose spatial integral gives the integer-valued baryon number, is given 
by 

24n^B = -EijktT {R^RjRk) ■ (3.12) 

From the mathematical point of view these field configurations represent 3-dimensional topo- 
logical solitons. 

The model is not integrable and thus, with few exceptions, explicit solutions are hard 
to obtain. One way to overcome this problem is by expressing the Skyrme field in terms of 
harmonic maps of 5*^ to CP^~^. In particular, the idea of the rational map ansatz, introduced 
in [13] is to separate the radial and angular dependence of the Skyrme fields. Its SU{N) 
generalization introduced in [16], expresses the Skyrme field in terms of a profile function g{r) 
and a N X N Hermitian projector P that depends only on the angular variable as 



[/(r,^,0)=e2*3M(^~i/^). 



(3.13) 



The matrix P is a harmonic map from 5*^ into CP'^ ^ . Hence it is convenient to map the sphere 
onto the complex plane via a stereographic projection in terms of the complex coordinate 



z = e^"^ ia.n{9 / 2) and its complex conjugate. In fact P = {V ® V^^)/|V"p is given in terms of 
a A^-component complex vector dependent on z and z. For fl3.13p to be well-defined at the 
origin, the radial profile function has to satisfy (7(0) = vr while as r — * cxo, it is required that 
the Skyrme field has \iv[Vj.^aog{f) = 0. 

In [H] it was shown that ansatz (13.131) describes field configurations for the SU{2) model 
that are close to being solutions of the model. In particular, although the corresponding 
energies are slightly higher that the energies of the exact solutions (obtained numerically) 
the symmetries of the baryon and energy densities are the same. Also in [16] , it was shown 
that when harmonic maps from S"^ to CP^^^ of the form (I3.13P are considered, low-energy 
configurations of the SU{N) Skyrme model can be derived. These configurations are more 
symmetrical than the SU{2) ones but have higher energies. However, exact solutions with 
spherically symmetry may also be obtained from (I3.13p . They correspond to solutions of the 
SU{2) and SU{3) Skyrme model with baryon number B = 1 and B = (topologically trivial 
solution), respectively. 

Using (13.130 the energy of the Skyrme model simplifies to 

^kin = ^ I g\ANr^ + 2Msin^g)dr (3.14) 

611 J 

Epot = ^j(ANr^gl + 2N{gl + l)sin^g + l^-^\dr (3.15) 



where 



An = ^(N-l) (3.16) 

AT = — [ dzdztr (\P,\'^) (3.17) 

2n J 

I = -^ f dzdz{l + \z\yti{[F,,P^]^). (3.18) 

Note that the integrals Af and X are independent of r. In particular. A/" corresponds to the 
energy of the 2-dimensional CP^~^ sigma model and is equal to the degree of the highest- 
order polynomial in z among the components of V^ = R{z) (when being holomorphic) after all 
their common factors have been cancelled out. 
Finally, the baryon number for this ansatz is 

B = ^ fdzdzti{F[F„F,]) (3.19) 

10 



which is the topological charge of the 2- dimensional CP^^^ sigma model. 

In general, for V holomorphic, it was proved in [16] that B = JV = n and 1 = n^ where 
n = deg(l^). However, for the special SU{3) non-topological solution the vector ^ is a function 
of z and z (i.e. non-holomorphic) . In this case, the baryon density is identically zero {B = 0) 
and the solution describes a bound state of two skyrmions and two antiskyrmions and as 
such is unstable as it corresponds to a saddle point of the energy. Let us emphasize that the 
field configuration is a genuine (non-trivial) solution of the SU{3) Skyrme model; while the 
corresponding parameters are A]\f = | and A/" = X = 4. 

The Bogomolnyi-type argument for the Skyrme model gives the following lower energy 
bound 




2X{gr + l) sin^g 

-2 (^/A^ + 2Af^ sin5f(9^(cos5f)|rfr 
> ^(2Af+^/A^y (3.20) 

due to the boundary conditions of the profile function. 

Baby skyrmions are topological solitons of the field theory which resembles the Skyrme 
model. Thus, in the following sections, we apply the techniques developed for discretizing the 
baby Skyrme model in the 3-dimensional Skyrme model. 

3.1 Discrete Skyrmions 

In what follows we present a lattice version of the Skyrme model by using the lower bound 
of the energy. Thus, following |3] we start with the same function smgdr{cosg) as appears in 
(I3.20p and reconstruct the inequality 

\/ A^I sing A {cos g) = —Dngn (3-21) 

where Dn — > \/ A^rg^ and F^ — > — sin^gf in the continuum limit h -^ ^. The formula 
A(cos5') = — 2//isin (^^^y^) sin (^^^*j-^) suggests the choices 

n nrt u\'^ ■ [9+- 9 

J^n = y A]\f[nn) — sm 



h \ 2 
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(3.22) 



Also, the origin must be treated in a special way since fl3.22p are undefined at n = 0. One 
possibility is to arrange it so Dq + Fq = implying that 



Dn = -Fn 



Vh 



■ \/9{.h,t) -vr sin g{h, t) 



(3.23) 



which follows for the discretization of the term: \/ AjqT [A((7 sin^ g) — gd>.{siT? g)\ . Note that, 
a direct discretization of (13.211) when ra = and assuming that Dq = —Fq is not possible since 
the terms are identically equal to zero due to the boundary condition g{nh,t)\n=o = vr. Thus, 
one possibility to overcome the problem is to discretize its counterpart term dr{g sin^ g) — 
gdri^sin^ g). This way the non-trivial ansatz (13.231) is obtained. 

Then, the potential energy of the lattice Skyrme model is given by 
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Finally, the discrete version of the kinetic energy is 
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(3.25) 



The corresponding Euler-Lagrange equations obtained from the Lagrangian L = ii^km — -E'pot 
where the potential and kinetic energy are, respectively, given by (13.241) and (I3.25p . read 
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(3.27) 



In the next section, we will show that stable discrete skyrmion solutions of the aforementioned 
equations can be obtained numerically. In particular, the analogues of the SU{2) B = 1 and 
SU{3) B = skyrmion configurations are constructed when the corresponding parameters are 

Aj^ = 1, J\f = 1 = 1 and A^ = 4/3, A/" = J = 4, respectively. 

3.2 Numerical Solutions 

The numerical existence and stability computations have been repeated similarly to the 2- 
dimensional case. The principal finding in this setting, as well, is that the skyrmion structures 
are found to be linearly stable, both in the SU(2) and in the SU{3) cases. In particular, 
examples of the profiles of the obtained structures and their linear stability are illustrated 
in Figure H] for values oi h = 0.4 (left panels) and h = 1 (right panels) for the two different 
parameter sets [top and bottom]. Notice the absence of any real eigenvalues showcasing the 
stabihty of the waves. 

Using numerical continuation, with the spacing h as the relevant parameter, reveals that 
the structures remain stable as h is modified. The monotonically decreasing trends of the 
relevant principal imaginary eigenvalues with h are shown in Figure \5\, together with the po- 
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Figure 4: Same as Figure [H but for the 3-diniensional case for a) the one SU{2) skyrmion and 
b) for the topologically trivial SU{3) skyrmion-antiskyrmion configuration. In each case the 
spatial profile of the skyrmion and the eigenvalues of its linearization are shown for h = 0.4 
(left panels) and h = 1 (right panels). 



14 



,,--13 




^- OA 



^"0.4 




Figure 5: Same as the panels of Figure[2]but for the 3- dimensional case when An =T = M = 1 
(top row) and An = 4/3 and X = M = ^ (bottom row). Notice that the right panels show 
the three lowest eigenvalues of the linear stability in each case. 

tential energy dependence on the spacing. The latter quantity has a non-trivial non-monotonic 
dependence with a minimum around a given spacing (which is dependent on the parameters 
Ani M and X, but is both cases occurring near h = 0.5). 

Finally, the numerical bifurcation analysis results were tested against the direct numerical 
integration of equations (13.271) . These simulations are once again particularly relevant in this 
context as the dynamical evolution equation contain terms (such as the second term in the 
left hand side of (I3.27P ). which are not accounted for at the linear stability level [arising at 
0(e^)]; hence, it is important to check whether linearly stable solutions may be destabilized by 
such higher order effects. Our results for these simulations are summarized in Figure El where 
integration results are shown for times up to t = 150, with a strong perturbation to the exact 
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Figure 6: In a format similar to that of Figure [3l the figure presents direct numerical simula- 
tions of the cases with ^^v = X = A/" = 1 (top row) and An = 4/3 and X = A/" = 4 (bottom 
row). The left panel shows in both cases the evolution of a skyrmion solution with h = 1 in 
space and time, while the right shows the deviation 6g{n, t) = g{n, t) — g{n, 0) from the initial 
profile. 

solution (for h = 1) oi magnitude 10~^ being imposed as the initial condition. However, the 
relevant perturbation remains bounded for the duration of the simulation and even for times 
up to twice as large as the ones shown here. These results clearly indicate the robustness of 
the obtained solutions. 

4 Conclusions 

In this paper a novel discrete version of the Skyrme model in 2 + 1 and 3 + 1 dimensions is 
presented. It has been shown that both models admit discrete skyrmion solutions, similar 
to the continuum ones, which are well-behaved and remarkably stable, in the sense that 
they cannot be destroyed by small (or even not so small) perturbations. The discretization 
scheme is based on using polar coordinates and thus, the corresponding skyrmions are radially 
symmetric; and therefore, it is not possible to investigate their fully 2d or 3d dynamics. The 
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advantage of the version described in this paper, on the other hand, is that the lattice spacing 
can be relatively large, without compromising the stability of the solitons. 

It should be instructive to produce discretizations of the full 2- and 3-dimensional systems, 
that are not restricted to radially-symmetric configurations, and to investigate the construction 
of the corresponding skyrmions. That way, the full dynamics of the discrete skyrmions can be 
considered, and directly tested against numerical simulations of the continuum ones [TT t fTS l fTQ]. 
The details of this program constitute an interesting direction for future study. 
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